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THE EFFECT OF REMOTE ZONES ON THE ACCURACY 


OF EVALUATING THE MOLODENSKY INTEGRAL 


James J. Buglla 


An analytical method is presented for determining^ 
the error in computing the gravity anomaly from 
the Molodensky integral if geoidal undulation data 
in remote zones are neglected. The error is given 
in terms of the usual degree variance and a set of 
parameters which are functions only of the size of 
the area in which undulation data are taken. A 
numerical calculation using Kaula's degree variances 
indicates that it is necessary to integrate the 
Molodensky integral over approximately a hemisphere 
in order to reduce the effect of neglected data to 
the order of 1 milligal in the computation of 
gravity anomaly. 


INTRODUCTION 

The Molodensky integral, which relates the gravity anomaly at a specific 
point to a global distribution of geoiaal undulation measurements, has 
long been of considerable theoretical use to the geodetic eonununity. 

This equation has recently begun receiving renewed practical attention 
because direct geoidal measurements are now available, through altimetry 
dnta obtained from the GEOS-3 spacecraft, thus permitting the direct 
Cv-'lcuiation of the anomalous gravity field. 
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Theoretically, the use of the Molodensky equation to compute gravity 
anomalies requires a continuous, global distribution of geoidal undula- 
tion measuronents. This ideal situation is, of course, never recdized 
in practice. The integral must be replaced by svaa over finite areas, 
and interpolation and/or extrapolation of measured data are generally 
required to produce a sufficiently dense network of data. 

It is thus of current Interest to Inquire if one can, as in the case of 
the Stokes equation, replace the global integration by an Integration 
over a small spherical cap centered at the computation point. Two related 
questions are thereby raised: What is the expected error incurred in 
the computation of the gravity anomaly by neglecting undulation measure- 
ments in the remote zone outside a spherical cap of specified size; and, 
how large should the spherical cap be in order to contain the error within 
a specified accuracy? 

A study which addresses these questions has been made and the results are 
the subject of this paper. The approach used is similar to that taken 
for addressing these questions for the related Stokes problem. An expression 
for the mean-square error incurred in neglecting data beyond a spherical 
cap of given size is derived in terms of a set of universal functions of 
only the cap size, and of the degree variances associated with the gravity 


field. 
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An illustration is provided indicating the mean-square error of omission 
for tvo different sets of degree variances, one using Kaula's rule of 
thumb, and the other for the GEM-U gravity model. It is shown that a 
cap size of approximately one hundred degrees in radius is required to 
reduce the mean-square of cmiission to about 1 mllligal. 

AhALYCTS 

The Molodensky integreJ. (Molodensky, et al., 1962; Pick et al, 1973) is 
generally witten as 

s p 

in which the reference figure to which all quantities are referred is an 
ellipsoid of revolution with equatorial radius, a, and flattening, f. 

Agp and are the gravity anomaly and geoidal undulation, respectively, 

at the point p, y and are mean values of the acceleration of gravity 

and radius of the reference ellipsoid, and do is the differential area 
element. N is the variable geoidal undulation associated with do and 
is measured along the normal to the reference ellipsoid, r^ is the chord 
distance between point p and do, and the subscript S on the integral 
denotes a global integration. 


k 


Equation (l) relates measured geoidal und\ilatlons about the reference 

ellipsoid, II and N , to gravity anomalies calculated vith respect to 

P 

the ellipsoid. Hence, all quantities in (l) should properly be computed 
on the ellipsoid. The chord distance, r^, euid the differential area, 
do, are in this case rather complex functions of position on the ellip- 
soidal surface. Fortunately, however, the Earth is very nearly spherical 
in shape with a flattening of approximately 1/300. Consequently, the 

expressions for r and do can be replaced by their spherical approxi- 
P 

mations vith an acceptably small loss in accuracy. Thus, if t is defined 
as the central angle between the point p and the area element do, then 


V 

r 2 R sin — 
p e 2 


Further, if p is taken as the pole of a spherical coordinate system, and 
a is the azimuth angle of do referred to any convenient meridian through 
P, the area element do can be written 


do R ^ sin Vd^do 
e 
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As a final preliminary, n n».*w function f(y) is defined by 


f(H') 



such that 



( 2 ) 


Equation (]) Is now written ns the sum of two parts - nn luicv'.ral over n 
spherical cap of anpiular radius centered at the point t , and h 

second inteRrnl eoveriru': the rest of the globe. 


'V 2-n 
n 



4'=H' a«0 

o 


f(T)(N - N ) slrtt' d'F di 
P 
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The second Integral of (3) gives the error that would be incurred in the 

computation of Lg by neglecting the undulation data in the region 
P 

beyond the spherical cap of radius We call this integral and 

refer to it as the "error of omission." It is the evaluation of c , and, 

P 

more importantly, the global RMS value of c , which is the topic of the 

P 

following analysis. 

Introducing a new function f(¥) with the follo%ring properties: 


f(t) = 0 for Y < V 

o 


f(’tf) = f(H') for t > Y , 

o 


the expression for becomes 


e 

P 



/ 


f('F)(N - N ) siny dH* da 
P 


(U) 


where now the integration covers the entire globe. 

Since f(4') is piecewise continuous, it can be expand.jd in a series 
of Legendro polynomials (Macmillan, 1958, p 386; Heiskanen and Moritz, 

1967 , p 28). 




00 
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f(’T) - ^ R^P^(co«t) (5) 

n«0 

where the constants are defined by [(2n ♦ l)/2]R^ to simplify the algebra 
later, and P^(cosV) are the Legendre polynomials with argument cosH*. The 
constants arc evaluated by utilizing the orthogonality relations for the 
Legendre polynomials (Heiskanen and Moritz, 1967, p> 30). 


2n»l 

2 


R 

n 


Tf 


/ 


4'=0 



f(4')P (co3V)sinY dy da 
n 


( 6 ) 


or 


R = 
n 


I 


y=o 


f(y)p (cosy)siny ly 

n 



f(y)p (cosy)siny dy 

n 


( 7 ) 


Substitute (5) into (U), interchange the order of summation and integration. 


and expand the integrals. 


e 

P 



n*0 




N P (coBr)slny dV da 
n 



P (co8Y)8iny dy da 


H can also be expanded in a series of surface harmonics 


n»2 


Where the N are given by equations similar to eq (6) of this text 
n 


n 2it 


2n+l 


n 


/ / 


(cosy)slny dt da 


y»o o»o 


Hence, the first Integral of (8) is simply 
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The second integral of (8) Is also easily evaluated by recadling that 

P (cosY) ■ 1, and by making use of the orthogonality relations 
o 



P (cosY)P (cosY)8inY df do ■ 0 , n m 

n m 

Uir 

“ 2n4l • " ■ “ 


and, henc'i, the second Integral is equal to Uw. 
Eq (6) is vritten in final form, then, as 



( 10 ) 


Eq (lO) gives the error in the calculation of the gravity anomaly Ag^ at 
a specific point P, caused by neglecting the undulation data beyond a 


spherical cap of anguleu* radius centered at P. Thus, is a 

function of the coordinates of p. 


While this is of considerable interest in its own right, a more useful 
parameter is one which has global applicability and is not identified with 
a specific point. 

One such parameter is the global RMS values of C , and is found by 

P 

squaring (lO) and evaluating the mean of the result of the entire globe. 
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By definition, M j^p^j value of integrated over the xinit 

sphere, divided by the curea of the unit sphere, or 



If we write N ■ x ^ 
p n 

n«2 


then e can be written 
P 


c 

P 


Sr ^ ■ "n’ "n 

n*2 


and thus (11) is 


M 




E 


00 

E 


(R - R )(R - R ) 
o n o m 


1 



N N do 
n m 


( 11 ) 


( 12 ) 


The bracketed term of (12) is the mean value of the product of two Laplace 
harmonics. This integral vanishes unless n ■ m. Heiskanen and Horitz 
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(196?) gives the following relations (pp. 2^6 and 2$9)t and show that the 
mean value of the product of two gravity anonaly harmonics can be written 
in terms of the degree variances - i.e. 

M lAg Ag I ■ 0 , n I* m 

I n ml 

■ , n ■ m (13) 


Also (p. 97 )« the general undulation harmonic can be written as a function 
of the gravity anomaly harmonic as 


d 

n 




(iM 


Using (13) and (1.1*), then, we can write 


M 


N 

n 



0 



Y2 



n m 


n ■ m 


12 


and (12) ccm be written in final form 






(15) 


This is the desired equation and expresses the mean square value of the 

error of omission in Ag , in terms of only the degree variances, c , 

P n 

and the terms, which are functions only of the cap size, 


NUMERICAL RESULTS 


Expressions for R can be found from a direct integration of eq (7) 
n 

for any y . R and R, are thus found to be 
o o 1 


"o ' ) 


R=R ~8(i-r) 

1 o 

Where F = sin(y /2). For larger n, the direct integration of (7) 

o 

becomes cumbersome. However, Meissl (1971) has derived a recursion formula 
for Moloden sky- type kernels. With a slight modification of notation this 


formula is 


13 


n 


2 ( P ^ (cost ) - P (cost )] 
n-2 o n o 

(n - 1/2) 8in(f /2) 
o 


♦ 2R - R ^ 
n-1 n-2 


vedues of R^ are presented in table I for n * 0 to 20, and for 

4*^ = 5° to l 80 ° in 5-degree increments. R^ through R^ are plotted 

against 4* in figure 1. All of the R tend towards zero as Y goes 

o no 

to l 80 ®, naturally, since the error of omission is zero IT one integrates 

eq (l) over the entire globe. All of the R^ display a rather heavily 

damped oscillatory character, with R crossing the zero-axis n times. 

n 

For moderate to large values of 4* , the magnitude of R decreases 

o n 

rapidly as n increases, indicating that the high frequency components 

of the gravity ''fs-ld contribute little to the error of omission for these 

cap sizes. For small values of 4' the contribution of the higher frequency 

o 

terms appears to be more significant. However, the degree variances them- 
selves drop off quite rapidly with increasing n. For example, using the 
well-known Kaula rule of thumb (Kaula, 1962), one can write 


c 


n 


96 (n - l)^ (2n + l) 

II 

n 


for the Earth's gravity field, and it can be seen that the c^ decreases 

approximately as 1/n. It has been found that, for 4* greater than 

o 

about 10°, seven or eight terms of (15) are sufficient for convergence to 
less than Ijf error in M Ie ^ !. 


lU 

For illustrative purposes, eq (15) was evaluated for two sets of degree 
variances — the Kaula approximation above, and a set of degree variances 
derived from the GEM-U gravity field (Lerch, et al., 1972). These results 


are displayed in fig. 2, where the quantity 


\/m ! 


I Gp j is plotted for a 


range of 4*^. The general trend of the curves is that expected, namely, 
that the error of omission drops off rather rapidly with increasing size 
of the spherical cap, asymptotically approaching zero error as goes 
to l80°. Both sets of degree variances yield quantitatively similar 
results, indicating that a cap size of the order of 90° - 100° is necessary 
to reduce the error of omission to the order of 1 milligal. 


CONCLUSIONS 


A rather simple formula has been derived which predicts the mean-square 
error incurred in the computation of gravity anomaly by the Molodensky 
formula when one uses only measured data in a small spherical cap surround- 
ing the computation point, and neglects data outside this cap. If one uses 
representative degree variances for the Earth it is found that integration 
of geoidal undulation data over approximately a hemisphere centered at the 
computation point is required to yield the gravity anomaly at the computa- 
tion point to an accuracy of 1 milligal. It is evident that rather extensive 
interpolation and/or extrapolation of measured data would i general be 
necessary to provide a sufficiently dense network of data for direct utiliza- 
tion of the Molodensky integral. 
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Table 1(b) P coefficients for ^- = 3 to l80 degrees 
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Table l(d) R coefficients for if; = 5 to 180 degrees 
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